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ABSTRACT 

The new interatomic potential g(r), providing the tunnelling of atoms, is sug- 
gested for calculating HFB phonon-roton spectrum E(p) of liquid 4 He. The 
tunnelling enables the calculations to reproduce the experimental dependence 
of shifts E(p) on pressure P induced by 3 He admixture in the constant total 
density and volume. 



Distortion of the excitation spectrum of liquid 4 He by admixture of 3 He was investi- 
gated experimentally in [1]. The authors of [1] appealed to the theoreticians to calculate 
the excitation spectrums in liquid-helium mixtures in order to explain their data. The 
first attempt [2] to explain the 4 He spectrum shift theoretically was based on the primary 
experimental results [3] , that were revised in [1] however. The recent theory [4] suggested 
a new version of interatomic potential and presented the calculation of 4 He spectrum 
shift caused by admixture of 3 He. Still the result [4] disagreed with the experiments [1], 
the explanation of that could be supposed in the invalidity of Hartree-Fock-Bogoliubov 
approximation (HFB) in [4]. Now we try to improve the situation using a more realistic 
potential than that in [4] in the framework of HFB approximation. We show that the 
experimental data can be explained by postulating that the realistic potential curve dif- 
fers from the usually assumed potential curve: the realistic potential has a local repulsive 
maximum in the attractive region of the usual curve. Moreover, the potential depends of 
the pressure (see Fig.l). 

In this paper, we present an interaction potential U, providing the tunnelling of atoms 
through the local maximum. This property is inherent for the atoms in liquids and fills 
the gap between discrete and continuous spectrums of ordinary two-atomic potential. 

Following [5], we assume that the polarization effects are not important in the liquid 
4 He and calculate the interaction potential U between two 4 He atoms with electrons in 
the ground state. We use the formula for the interaction energy of two atoms with charge 
distributions densities pi(r) and pi{v) 

u= f Pi(r)p 2 (r')drdr' ^ 
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where the charge densities are 

pi(r) = q5(r - n) - qF(r - n), 

(2) 

p 2 (r>) = q S(r' - r 2 ) - qF(r' - r 2 ). 

Here F(r — ri) and F(r' — r 2 ) are the electron probability densities at the points r and r', 
the atom nuclei being at the points r x and r 2 , and q = 2e. The interaction between atoms 
contains three terms: U nn , U nc , U cc , describing interactions between two nuclei, nucleus 
of one atom with the cloud of the another and between two clouds. 
The term U nc has the form 

2 f 5(r - ri)F(r' - r 2 )drdr' 2 r F(r' - r 2 )dr' 2 r F(r )dr 

Unc = ~q / 1 7; = -q / — | 7j — = -q / 75 r, (3) 

J \r — r'\ J |ri — r I J |R — r | 

where r = r' — r 2 and R = r x — r 2 is a 3D vector between nuclei. The electron probability 
density F(r) is determined by the electron wave function $ of the 4 He atom. Therefore 
we use the Hylleraas-type wave function with parameters 7^2 (see [6]) 

$ = e- s/2a [I+71-+72V! > s = r 3 + r 4 , 
\ a a A J 

t = r 3 - r 4 , it = |r 3 - r 4 | , 

where r 3 and r 4 are the electron coordinates of atoms a is the Bohr radius. In order to 
obtain the probability density (2) of the cloud F at the point r one can add the W\ and 
W2 probabilities to find one electron at the point r, and another at an arbitrary point 

qF(r) = ewi(r) + eu^r), 

wi(r) = J |$(r 3 = r,r 4 )| 2 dr 4 , w 2 (r) = J |$(r 3 ,r 4 = r)| 2 dr 3 , 
These expressions show, that the probability density F(r) can be written as 

F{r) = -±- e - 2r / a (a 4 r 4 + a 3 r 3 + a 2 r 2 + a x r + a + a^~) = -^e- 2r ^. (4) 
ira 6 r -na- 

Thus, the formula for U nr is 



q 2 f e~ 2ro ^{r Q )dr 

U nc = 3 / P5 ; • 5 

71CL J |R — r | 



The direct computations give 

U = const • e~ 2x ■ (fox 6 + (3 5 x 5 + fax A + (3 3 x 3 + 



1 \ R 

+(3 2 x 2 + (3 1 x + (3 + P-x- ) > x = —. (6) 



The parameters a in (|j) are the functions of parameters (3 in ®. We emphasize that the 
potential (BJ) decreases exponentially at the distances R^> a. 
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Fig. 1: Initial (1) and shifted (2) potential g(x). 

For the sake of simplicity, we limit ourself by considering the potential of the form 

g (x) = const ■ e~ 2x (Ax 6 + Bx 5 + C + D^j , (7) 

instead of U(x); g(x) has the local maximum and depends on the minimum number of 
parameters (see Fig. 1). 

It was shown in [1], that the pressure P increase in 4 He due to the 3 He admixture (in 
the same volume V and for the same total atomic density) leads to the deformation of 
phonon-roton 4 He excitation spectrum E(p). Our purpose is to reproduce the deformation 
A(E) = E(p)\p — E(p) of spectrum curves caused by the pressure P. The HFB formula 
for the spectrum of noncondensate excitation is 



E(p) 



p 4 , p 2 . 

—pg{p), 



Am 2 m 



where m is the mass of atom, p is the condensate density and g is the Fourier transform 
of g. The increase of P leads to decrease of p and to the deformation of the electron 
clouds of the atoms, which keep, however, spherically symmetric. This is evident that the 
cloud deformation results in the displacement of the potential minimum go "down" and 
"to the left" in the x, g(x) plane. These curves 1,2 are shown in Fig. 1 and corresponds 
the values 

A 1 = -0.5, B 1 = 1.4, Ci = -2.5, D 1 = 0.2, 



A, 



0.2. 



-0.59, B 2 = 1.8, C 2 = -3.8, D 2 

The choice of parameters A,B,C and D is determined by the conditions for the potential: 
to be positive (repulsive) for x <^ 1 (D > 0), 
to be negative (attractive) for x 3> 1 (A < 0), 
to have local maximum. 
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Fig. 2: Shift of E(y) with the change of parameters: A = -0.5, B = 1.4, C = -2.5, £> = 
0.2, p = 0.03 (solid line); A = -0.59, 5 = 2.3, C = -19, L> = 0.2, p = 0.027 (dashed 
line); A = -0.61, B = 4.5, C = -30, D = 0.2, p = 0.025 (dotted line). 



The local maximum of g(x) occurs between coordinates of go and g' . This maximum 
is introduced in order to describe the tunnelling of atoms, providing the agreement of the 
present calculation with the experimental results for A(E) in [1]. This tunnelling was not 
taken into account in [4], and this is the reason for a failure of [4] to explain [1]. 

In order to calculate the spectrum, we need the Fourier transform of g(x) 



/ \ f Sm X U I \ 2 7 

g(y) = —const ■ a / g{x)x ax = 

o " 



—const- a 3 



AT (8) (4 + y 2 ) 4 sin [8 arctan(y/2)] + BT(7) (4 + y 2 ) 7/2 sin [7 arctan(y/2)] + 
AC D 1 

+ Ta , + ta \ 2\ ' where y = a P- 

(4 + y 2 ) (4 + y 2 )_ 

We let go and g' Q denote the minima of the potential U(x) (see Fig. 1), and consider the 
behavior of U(x). The curves 1 and 2 in Fig. 1 demonstrate the general situation. 

In order to obtain the experimental scale of spectrum shifts we use the Fourier images 
of potentials similar to Fig. 1, but different from them in scale. Energy excitation spectra 
for such potentials are shown in Fig. 2. The shifted curves (dashed, dotted) in Fig. 2 
correspond the decrease of condensate density p and increase of pressure P, relative to 
the initial state (solid curve). The shifts of the curves in Fig. 2 qualitatively correspond 
the deformations of g (x) shown in Fig. 1 . 

The two types of shifts of initial curve are possible. If the potential has a local 
maximum, then there are no intersections of E(p). We emphasize that if the potential 
has not a local maximum, then we obtain the intersections of the shifted spectrum curves 
as shown in [5]. 
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The local maximum implies the tunnelling of atoms and we conclude that this prop- 
erty is inherent for quantum liquid. It means, that in the case of ordinary Lennard- Jones 
potential there are discrete and continuous parts of spectrum only. In our case the tun- 
nelling of atoms is possible in the range of energies intermediate between continuous and 
discrete spectrum. This tunnelling enables us reproduce the experiment sign of A(E). 
It describes the atoms, penetrating in liquids through the energy gap corresponding the 
regimes of "nearest atoms ordering" (in solids) and "atomic distant disorder" (in gas). 
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